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In  this  paper  Jahn-Teller  vibronic  coupling  calculations  are  presented 
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for  the  general  case  in  symmetry.  The  Hamilton i^^|^loyed  is 

complete  up  to  terms  quadratic  in  vibrational  coordinates  qf  and  q which  . 
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It  is  concluded  that,  in  agreement  with  previously  reported  spectroscopic 
data  for  ReFg  and  IrFg,  quadratic  terms  in  the  vibronic  interaction  are 
essential  (both  qualitatively  and  quantitatively)  to  the  complete  under- 
standing of  intrastate  vibronic  coupling  in  transition  metal  hexafluorides. 
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I.  INTRODUCTION 


1 -T 

Recent  experimental  studies  of  transition  metal  hexafluorides  have 
indicated  that  the  spectroscopically  observed  dynamic  Jahn-Teller  (JT)  effects 
are  not  adequately  explained  by  the  usual  linear  JT  theory^;  quadratic  terms, 
at  least,  are  required  in  the  vibronic  Hamiltonian.  In  order  to  treat  hexa- 
fluoride data  quantitatively,  it  would  be  desirable  to  solve  the  problem  of 

two  JT-active  vibrations  (v2  (e  ) and  (t,_))  vibronically  coupled  to  a 

9 ^9  oq 

electronic  state  for  a Hamiltonian  complete  up  to  quadratic  terms  x (e^  + 
^2g^^QUAD'  Unfortunately,  however,  the  requisite  numerical  methods  are  not 
feasible  given  currently  available  computers.  In  fact,  even  the  [Pg^  x t2g]g^j^Q 
problem  alone  presents  an  arduous  task.  The  [Pg^  x coupling  problem 

is,  however,  readily  soluble.  Its  solution  is  presented  here  and  is  found 
to  generate  a qualitative  insight  into  the  nature  of  the  quadratic  portion  of 
the  JT  interaction. 

Previous  work  on  [Pg^  x Sglgy^g  (or  the  very  closely  related  problem: 

4 5 

^^9  ^ ®g^0UAD^  focused  on  the  weak  coupling  limit  or  on  the  case  for 

which  linear  coupling  is  large.®  In  the  present  work  matrix  elements  are 
given  which  allow  full  secular  matrices  for  [Pg^  x Sgjgy^g  to  be  set  up  and 
diagonalized  numerically.  Given  a sufficiently  large  enough  truncated  secular 
matrix,  accurate  eigenvalues  and  eigenvectors  can  be  obtained  for  any  values 
of  the  linear  and  quadratic  parameters.  These  secular  matrices  are  also  appli- 
cable to  the  [Eg  x ®g]gu/^[)  situation. 


iT^' 
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II.  THEORY 


The  first  step  toward  obtaining  the  necessary  matrix  elements  for 
^ ®g^QUAD  vibronic  Hamiltonian.  Englman^  has  discussed 

in  detail  the  group  theoretical  techniques  which  allow  a parametric  Hamil- 
tonian to  be  obtained;  it  is  given  here  with  only  minor  modifications: 


^ y are  the  mass-weighted  normal  coordinates  of  the  eg  vibration, 

^ p ^ are  the  conjugate  mass-weighted  momenta,  is  the  vibra- 
tional force  constant,  I^  is  the  4 x 4 Identity  matrix,  (t^,  C^,  C^)  are  the 
linear  and  quadratic  coupling  parameters,  and  ( ^ , are  Dirac 

matrices.^  In  the  following  discussion  we  will  also  make  use  of  o^,  oj^, 
which  are  additional  Dirac  matrices  as  given  In  reference  4. 

Choice  of  vibronic  basis  set  is  important  since  many  matrix  elements  need 
to  be  calculated.  In  any  case,  the  basis  will  Involve  two-dimensional  harmonic 


oscillator  functions,  {X(n2,i2)l  and  vectors  which  represent  the  four  elec-  | 

tronic  components  {^(03,  03)}.  The  {V(p3,a3)}  are  eigenvectors  of  both  and  1 

J 

^ with  eigenvalues  P3  » ± 1 and  03  = ± 1.  The  question  of  how  to  combine  ij 

these  factors  into  a vibronic  basis  set  that  is  most  convenient  is  now  addressed.  I 

One  approach, which  will  only  be  outlined,  involves  utilizing  molecular  point 
group  theory.  The  vibrational  factors,  however,  do  not  transform  as  standard 

j 

irreducible  representations  of  Ojj  but  appropriate  linear  combinations  of  these 

functions  can  be  generated  which  do  so  transform.  The  symmetrized  vibrational 

$ 

factors  and  the  electronic  factors,  which  transform  as  fg^,  can  be  combined  i 

using  Clebsh-Gordon  coefficients  to  give  vibronic  functions  which  transform  | 

as  Tc  , r,  , or  . Separate  secular  matrices  can  then  be  formed  in  the 

bg  /g  og  ^ 

usual  fashion.  This  method  will  not  be  used,  however,  because  the  calcula-  | 

1 

tion  of  matrix  elements  within  this  basis  set  is  more  cumbersome  than  need  | 

be.  I 

An  alternate  approach  is  to  use  the  same  vibronic  basis  that  is  used  in  | 

the  linear  problem.  These  are  written  as  follows:  I 

b {n2,  P3»  03;  J2)  = X(n2,t2)  V (p3,  03)  (6)  I 

in  which 

J2  = ^2  ^ ” ±1/2,  ±3/2,...  j 

n2-0,l,2,...it2~“^2*~^2'^^*****^2  I 

P3  = ±1 , O3  = ±1 . I 

In  linear  coupling,  J2  is  a good  quantum  number  and  separate  secular 

matrices  are  generated  for  each  value  of  J2.  Upon  introduction  of  quadratic  \ 

terms  (/-^  '(Cg)  in  particular),  J2  ceases  to  be  a good  quantum  number.  One  i 

can  determine,  by  examining  a general  matrix  element,  that  [ZJj  '’’oU  3}  remains  j 

a good  quantum  number.  Therefore,  the  following  three  sets  of  basis  vectors,  j 

labeled  by  J2,  are  coupled  by  the  quadratic  terms:  [ 


f 
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{±3/2,  ±9/2,  ±15/2...)^° 
{...,  -5/2,  +1/2,  +7/2, 
{...,  -7/2,  -1/2,  +5/2, 


It  is  of  interest  to  correlate  the  {2  J2  mod  3}  quantum  number  with  irreduci- 
ble representations  of  0^.  It  can  be  shown  that  {0  mod  3}  correlates  with 
r,  and  r,  while  both  (1  mod  3)  and  {2  mod  3}  correlate  with  Tq^.  Thus,  one  dis- 
advantage  of  this  basis  set  is  that  the  (0  mod  3}  block  is  not  factored  as 
much  as  possible  (separate  and  sets);  nonetheless,  with  modern  com- 
puters and  diagonal ization  routines,  this  is  not  a serious  drawback. 

The  main  task  remaining  is  to  find  integrals  of  powers  of  vibrational 
coordinates  over  harmonic  oscillator  functions;  these  are  often  referred  to 
as  the  primitive  matrix  elements.  General  formulae  exist^  which  allow  these 
matrix  elements  to  be  evaluated  simply. 

Primitive  matrix  elements  relevant  to  the  case  at  hand  are: 


/ X*(n2. 
/ X*(n2, 
/ XMnj, 
/ X*(n2, 
/ X*(n2. 


12)9+  ^('^2  + 1 » ^2  + = 

11  - » 1/2 
■ 1 (n2  + *.2  2)  1 , 

(7) 

12)^1  X(n2  +2,  e2  + 2)dT  = 

+ 

1 i )l/2 

^ j (n2  + ^2  + 2)(n2  + t2  + 4)J  , 

(8) 

ii2)q^  X(n2,  ^2  t ^ 

1 4(12  + t2)(n2  + il2  2)|  ^/^, 

(9) 

t2)9+9.X(n2  + 2,  «,2)‘^'c  = ^ 

{ (n2  - 1I2  + 2)(n2  + ^2+  2)j  1/^, 

(10) 

t2)q+9.X(n2.  t2)dT  = ^ {2n2  + 2}  , 

(11) 

in  which 


% = %t  \ 
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The  desired  vibronic  matrix  elements  are  then  readily  found: 


/b  (nj,  i2’,  P3.  03*t  ^("2  + 1.  t2  + P3;  -P3,  03;  Jajdt 


ii  oj 


{D2(n2  + P3e2  + 2)} 


fb 


(f'2»  ^2*  P3»  P3J  ^z)S^  ^^^(aT^)b(n2»  iz\  P3,  03;  J2)dT  = — + ]}  ^ (13) 

t»  U) 


/b  (n2,  ^2;  P3»  03;  J2  (aig)b(n2  + 2,  t2*  P3»  P3i  " 


tlu) 


QataiJ  V2 

--^--^-3-  {(n2  + t2  + 2)(n2  - ^2  + 2)}  , 


fb 


(n2.  ^2;  P3«  f’3*  J2)  ^^^^eq)b(n2,  *'2"2p3‘*  ■P3»  <^31  J2*3p3)dT  = 


% 


,1/2 


Q2[egJf4(n2  + P3*-2)(n2  ■ P3t2  + 2))  , 


/b  (n2,  ta;  P3.  03;  Ja)^  (e  )b(n2  + 2.  i2-2P3;  -03.  03;  J2-3p3)dT  = 

^ 0) 

c 


Q2[eg]  {(n2  - P3t2  ■*■  2)(n2  - P3II2  ■*■  ^)1 


1/2 


in  which  the  dimensionless  coupling  parameters  are  defined  as:  D2  ° 


^ ’ ^2[eg]  * ^ 


(12) 


(14) 


(15) 


(16) 
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Note  that  03  does  not  appear  on  the  right  hand  side  of  any  of  these  equations 
and  thus  two  Identical  secular  matrices  are  generated,  one  for  each  value  (+1) 
of  the  quantum  number  03.  It  will  prove  to  be  Interesting,  for  later  consider- 
ations (see  Section  III),  to  use  perturbation  theory  to  find  the  splitting  of 

A /pN 

the  r\2  = 1 level  due  to  unrealistic  case  of  D,  = 0 (no  linear 

coupling).  The  secular  determinant  which  results  from  a first  order  degenerate 
perturbation  theory  treatment  of  the  02  = 1 level,  employing  equation  15,  is 


(2  - X)  4 QaCeg] 
4Q2[eg]  (2  - X) 


(17) 


with  solutions 

X * 2 + 4 QaCCg]. 

These  quantities  are  measured  in  units  of  vibrational  energy,  It  Is  possi- 

ble  to  obtain  a first  order  expression  for  .i^'^^(aig)  also;  however,  this  is 
not  necessary  since  an  exact  expression  easily  obtains: 

= (n^  + 1)  >/l  + QaUig]  = (n^  + Dd  + 


X 


+...)  . (18) 
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III.  DISCUSSION 

Examples  of  the  secular  matrix  calculation  of  energy  levels  of  [Pg  x 
are  given  in  Figures  1 and  2.  The  secular  matrix  used  in  these  calculations 
was  truncated  after  r\2  = 10;  a basis  of  this  size  results  in  44  x 44  matrices 
for  both  the  {0  mod  3}  and  (1  mod  3}  blocks  (the  {2  mod  3}  block  gives  the  same 
eigenvalues  as  {1  mod  3}).  Under  these  circumstances  it  is  required  that  D2 
be  less  than  1 for  an  accurate  description  of  the  levels.  The  upper  limit  of 
QzE^q]  was  not  numerically  determined,  but  is  probably  about  0.5.  The  behavior 
of  the  levels  at  these  large  values  of  D2  and  Q2  [e^]  is  quite  complicated 
due  to  specific  level  repulsions. 

D2  may  be  taken  as  positive  since  the  relevant  matrix  elements  are  all 
off-diagonal  (Eqn.  12).  The  situation  with  respect  to  Q2  [Cg]  is  more  compli- 
cated even  though  it  too  has  all  off-diagonal  matrix  elements  in  the  employed 

A 

linear  JT  basis  (Eqn.  15,  16).  ^ (e  ) will  have  diagonal  matrix  elements 

in  a symmetry  adapted  basis  set.  The  eigenvalues  of  (e^)  are  inde- 

pendent of  the  sign  of  Q2  C®g]»  however,  the  concomitant  changes  in  the  eigen- 
vectors upon  changing  the  sign  of  Q2  [Cg]  may  alter  their  symmetry  transforma- 
tion properties.  For  Q2  [Og]  -*■  - Q2  [Cg].  Py  and  Py  Pe- 

Several  general  observations  can  be  made  based  on  these  calculations: 

i)  For  the  case  D2  = 0 (Figure  1),  the  perturbation  expression  Eq.  17  is 

found  to  be  useful  over  a large  range  of  QaCfig]  values  (Q2[eg]  <0.1)  by  direct 

comparison  with  results  of  the  truncated  secular  equation  calculation.  This 

^ (21 

situation  arises  because  '(eg)  couples  02  and  n2  + 2 levels,  but  not  n^ 
and  n2  + 1,  as  in  the  linear  case. 

ii)  In  the  regime  (D2  <0.2,  QaC^gl  ^0.05),  for  the  02  = 1 levels,  effects 
of  linear  and  quadratic  terms  are  found  to  be  approximately  independent;  that  is, 
the  (ny  = 1,  Jy  = 1/2)  level  is  not  shifted  appreciably  as  Q^iCCg]  increases,  center 
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of-gravity  of  the  (112  = 1,  J2  = 3/2)  levels  is  preserved,  and  splitting  of 
(n2  = 1,  J2  = 3/2)  is  approximately  the  same  as  in  the  D2  = 0 case. 

iii)  Behavior  of  the  02  = 2 levels  are,  however,  qualitatively  differ- 
ent from  those  of  02  = 1 (see  Figures  1 and  2);  one  difference  is  in  the 
behavior  of  the  (02  = 2,  J = 3/2)  levels  (Figure  1).  For  D2  = 0,  these 


levels  do  not  split  under  c^^^^(eg),  whereas  D2  / 0 allows  admixture  of 


n2  = 1,3,...  levels,, thus  causing  (n2  = 2,  J2  = 3/2)  to  split.  It  should  be 

-A 

noted  that  existing  perturbative  treatments  of  ^ are  not  able  to  generate 


this  splitting,^  although  a more  complete  perturbation  calculation  would,  of 


course,  reproduce  these  results. 


Based  on  the  form  of  the  perturbation  energy  expressions  for  the  linear  and 

A 

quadratic  JT  terms,  it  can  be  seen  that  certain  n2  levels  are  split  by 

in  the  first  order  of  perturbation  theory  Eq.  17)  while  the  same  02  levels 

8 


are  split  by  ^ only  in  the  second  order."'  It  is  thus  possible  for 
quadratic  JT  terms  to  be  more  effective  at  splitting  vibronic  levels  than 
are  linear  JT  terms.  Examination  of  the  perturbation  expressions  for  the 
n2  = 1 levels  verifies  this  idea: 


A (3J  = 2 1 2(0^)^ 


X (c^^^^eg))  = 2 + 4Q2[eg]. 


(19) 


is  used  for  comparison  with  Q2[®g]  since  it  is  proportional  to  just  as 


Q2[eg]  is  proportional  to  C^..  Thus,  for  D2  = Q2t®a3  " ^9S.  19  show  that 


the  splitting  induced  by  Is  an  order  of  magnitude  m^vo  than  that 


induced  by  Similar  results  obtain  for  <^^^^(aig)  (see  Eq.  18). 


■ 1 K 
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1 -3 

Available  hexafluoride  data  show  that  the  splitting  induced  by  the  linear 
terms  is  of  the  same  order  of  magnitude  as  the  quadratical ly-induced  splittings 

A 

and  shifts.  Although  these  data  indicate  the  expansion  of^^  in  powers  of  q. 
and  q is  not  converging  as  rapidly  as  might  be  desired,  the  entire  effect 

^ A 

cannot  be  attributed  to  non-linearity  of  . Greater  effectiveness  of 

) in  causing  shifts  and  splittings  is  also  an 

important  consideration. 
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IV.  CONCLUSION 

An  accurate  numerical  method  is  given  for  treating  vibronic  coupling  of 

fo.,  (OJ)  electronic  state  to  an  e vibration  for  a vibronic  Hamiltonian  with 
09  n 9 

linear  and  quadratic  terms.  Examples  of  the  calculation  are  given  and  general 
comments  on  the  behavior  of  the  levels  under  the  Influence  of  quadratic  terms 
are  made.  It  Is  also  found  that  the  quadratic  coupling  term  is  a more  effective 
perturbation  than  the  linear  coupling  term  with  respect  to  the  observed  spectro- 
scopic splitting. 
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Figure  1 . 

Quadratic  Jahn-Teller  calculation  for  coupling  of  a Tg  electronic  state  with 
a V2(®g)  vibration,  (Tg  x The  fixed  parameter  values  are  D2  = 0 and 

V2  = 670  cm"\  D2  is  the  linear  Jahn-Teller  coupling  parameter  and  vf  is  the 
unperturbed  \>2  harmonic  oscillator  frequency.  Q2  [e^]  is  the  quadratic 
coupling  parameter.  Note  the  large  range  over  which  the  splitting  is  linear 


same  secular  matrix  ( the  {0  mod  3}  block,  see  text).  The  symmetry  labels  of 


the  eigenvectors  are  generated  from  their  transformation  properties.  Changing 
the  sign  of  Q2  t^g]  will  not  alter  the  energies,  but  will  interchange  Fg  and 
symmetry  labels. 
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I 
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Figure  2. 


Similar  to  Figure  1 but  D2  has  been  set  to  0.125.  Note  that  the  n2  = 1 for 
this  value  of  D2  levels  behave  similarly  to  those  for  which  D2  = 0 in  that 
one  pair  of  levels  splits  symmetrically  and  the  other  level  remains  largely 
uneffected  as  Q2  [e^]  varies.  This  situation  does  not  hold,  however,  for  the 
n2  = 2 levels.  The  J2  quantum  number,  which  is  good  at  Q2^®g^  ” 


indicated. 
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